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The two dimensional Hubbard model in the presence of diagonal and off-diagonal disorder is studied 
at half filling with a finite temperature quantum Monte Carlo method. Magnetic correlations as well 
as the electronic compressibility are calculated to determine the behavior of local magnetic moments, 
the stability of antiferromagnetic long range order (AFLRO), and properties of the disordered phase. 
The existence of random potentials (diagonal or "site" disorder) leads to a suppression of local 
magnetic moments which eventually destroys AFLRO. Randomness in the hopping elements (off- 
diagonal disorder), on the other hand, does not significantly reduce the density of local magnetic 
moments. For this type of disorder, at half-filling, there is no "sign-problem" in the simulations 
as long as the hopping is restricted between neighbor sites on a bipartite lattice. This allows the 
study of sufficiently large lattices and low temperatures to perform a finite-size scaling analysis. 
For off-diagonal disorder AFLRO is eventually destroyed when the fluctuations of antiferromagnetic 
exchange couplings exceed a critical value. The disordered phase close to the transition appears to 
be incompressible and shows an increase of the uniform susceptibility at low temperatures. 



I. INTRODUCTION 

Electrons in crystals are scattered both by their mutual 
interaction and by static disorder potentials. These pro- 
cesses typically lead to quite different or even competing 
effects. For example, on a bipartite lattice close to half- 
filling the strongly screened Coulomb interaction between 
electrons can generate antiferromagnetic long range or- 
der (AFLRO) while disorder tends to destroy such cor- 
relations. The simultaneous presence of interaction and 
disorder cannot in general be considered as a simple su- 
perposition of both contributions but new many-body 
phenomena may emerge. This has been found for in- 
stance in the study of the metal-insulator transitions in 
doped semiconductors [jjj or in the stability of AFLRO 
against disorder within a dynamical mean- field theory §. 
In spite of the great progress that has been achieved in 
understanding interacting as well as disordered systems 
in recent years, there is still no controlled and at the same 
time tractable theoretical method to describe their com- 
bined effects, in particular when the interactions and/or 
disorder cannot be considered as small. 

It is the purpose of the present paper to provide 
approximation-free results for a very simple microscopic 
model that incorporates electron interactions as well 
as disorder, namely the disordered Hubbard model, or 
"Anderson-Hubbard model" in two dimensions, D = 2. 
The Hamiltonian of the model reads in the usual nota- 
tion: 



H= t ij 6 L 5 j (T + ^(ei-^)7ii (T 

<ij>,<r iff 

+f7^(n 1T -l/2)(n u -l/2). 



(1) 



i and j are lattice vectors. The distributions of the ran- 



dom hopping elements fy and random local potentials ei 
will be specified later on. Note that hopping processes 
are restricted between nearest neighbors on a square lat- 
tice, hence there is no magnetic frustration. 

Model (|l|) has been investigated within various ap- 
proaches: the formation of localized magnetic moments 
has been studied in D = 3 for a very broad distribution 
of tij within an unrestricted Hartree-Fock approximation 
||. The case of random potentials has been treated by 
a real-space renormalization group method in D = 1, 3 
Q and D — 2 ||. In D = 2 this treatment provides 
a transition from a Mott to an Anderson insulator with 
no metallic phase. However, in these investigations the 
formation of AFLRO which will set in, at least in the 
unfrustrated case close to half-filling, is not taken into 
account. The (in)stability of AFLRO with respect to 
diagonal disorder and two types of metal-insulator tran- 
sition were examined [0 in a dynamical mean-field theory 
which becomes exact in the limit of D — > oo Q . Diagonal 
disorder has also been studied in D = 3 by Hartree-Fock 
approximations. The strong-coupling limit of model 
(|l|) with diagonal disorder was studied Q using a slave- 
boson formulation of the corresponding t — J model. In 
the case of off-diagonal disorder at half-filling the model 
maps in this limit onto the spin-1/2 Heisenberg model 
with random (not-frustrated) exchange couplings. This 
model has also been investigated numerically in D = 2 
H . Finally, the one dimensional Hubbard model with ei- 
ther type of disorder was studied using quantum Monte 
Carlo (QMC) simulations fo|. 

Here we will concentrate on the effect of disorder on 
the magnetic correlations. We will address the following 
questions: 

• How arc short and long range magnetic correlations 
affected by the two different kinds of disorder? 
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• Is there a critical disorder strength where AFLRO 
ceases to exist? 

• How do the magnetic susceptibility and charge 
compressibility behave in the disordered state? 



II. COMPUTATIONAL METHOD 

We will study the Anderson-Hubbard model in D = 2 
using a finite temperature, grand canonical quantum 
Monte Carlo (QMC) method jnj which is stabilized at 
low temperatures by the use of orthogonalization tech- 
niques The algorithm is based on a functional- 
integral representation of the partition function by dis- 
cretizing the "imaginary-time" interval [0, 0\ where (3 is 
the inverse temperature. The interaction is decoupled by 
a two-valued Hubbard-Stratonovich transformation H] 
yielding a bilinear time-dependent fermionic action. For 
the positive U model, the D + 1 dimensional auxiliary 
field (sii = ±1 where i is the lattice and / the "time" in- 
dex) couples to the local magnetization (njf — rii\). The 
fermionic degrees of freedom can be integrated out an- 
alytically, and the partition function (as well as observ- 
ables) can be written as a sum over the configurations 
of the auxiliary field with a weight proportional to the 
product of two determinants, one for each spin species. 
The two determinants are not equal since Su couples with 
different sign to the two fermion species, and, in general, 
their product is not positive definite and thus cannot 
serve as a weight function in an importance sampling 
procedure. The formally exact treatment of this "minus- 
sign problem" can lead in some regimes of the model pa- 
rameters to very small signal-to-noise ratios in physical 
quantities that become in fact exponentially small with 
inverse temperatures and system size. 

In the case of a bipartite lattice, under the particle-hole 
transformation of one spin species: 



c iL — > (-1)%. 



(2) 



Hamiltonion (|l|) is mapped onto the negative U Hubbard 
model: 

#(E0— >#(-E0+5>i-/i)(l-2rm). (3) 



If the Hamiltonion is now spin-symmetric, i.e. e\ = fj, 
for all i, the two determinants for spin up and down 
are identical, since in the case of the negative U model 
the Hubbard-Stratonovich field couples to the charge 
(riii that is, with the same sign for the two fermion 

species. Hence the determinant product is positive (semi) 
definite and there is no "minus-sign problem". For this 
reason the local random potentials ei lead to a minus- 
sign problem even at half filling whereas the sign is al- 
ways positive at half filling for off-diagonal disorder as 
long as the hopping remains restricted between nearest 
neighbors. 



We will consider the case of static, uncorrelated dis- 
order in either the hopping elements iy or the on-site 
potentials e\. Therefore we have to average all quanti- 
ties over a sufficient number of disorder realizations and 
calculate the averaged expectation values: 



{{A)) h = II f / *y^(*y) 



(A) ({Us}) (4) 
» ({*}). (5) 



(A) denotes the thermal expectation value of the operator 
A for a given disorder configuration; (• • -)fc( g ) stands for 
the average over hopping (site) disorder. We will assume 
uniform distributions of either the hopping elements iy 
with an average value of t = —1 or the on-site potentials 
£i with an average of zero. The width of the distributions 
Ph(s) is denoted by A: 



Ph(t ii ) = ~9ij-\t li -t\) 



(6) 
(7) 



(We use the same symbol for the width A since each 
type of disorder is considered separately.) Restricting 
to the half-filled band case (/x = 0) the remaining three 
parameters are: interaction U, disorder strength A, and 
temperature T = 1//3. 

In the present study we will concentrate on the follow- 
ing observables: i) magnetic correlation functions: 



(8) 



Here rhj = J2 a £7 Ai<7 i s the local spin operator, and N 
is the total number of lattice sites. y/C(0, 0) measures 
the density of local magnetic moments and is equal to 
n — 2d with the electronic density n — 5Zj ff ((Aj(T)) an< ^ ^ nc 
density of doubly occupied sites d = J2j(( n i7 n il)) ■ (The 
indices h, s of the disorder averages are suppressed for 
convenience.) ii) magnetic structure factors, the fourier 
transformation of (7(1): 



S(q) = X;C(l)e iqI ) 



(9) 



(note that (3S(0, 0) is equal to the uniform spin suscepti- 
bility), iii) charge compressibility: 
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(10) 



with local charge operator nj =^2 n- ia . 

We found S(q) to be largest at the commensurate vec- 
tor q = (jr, it). At sufficiently low temperatures the 
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magnetic correlation length exceeds the system size and 
S(ir, 7r) saturates with (3 for a given system size. Using 
the saturated values we can extrapolate the behavior in 
the thermodynamic limit by a finite size scaling accord- 
ing to spin wave theory Jl5| which predicts in the case of 
AFLRO in the ground state with sublattice magnetiza- 
tion M : 



C{N x /2,N x /2) 
S(ir, 7r) 



OiN- 1 ) 



N 



T 

M 2 

-Y + oiN- 1 ), 



(11) 



where (N x /2,N x /2) is the maximal separation on a 
square lattice of linear size N x = yN with periodic 
boundary conditions. Thus, we have two independent 
quantities to extrapolate the value of the ground state 
order parameter. The finite size extrapolation is tech- 
nically only possible in the case of off-diagonal disorder 
where there is no "minus-sign problem" and hence suf- 
ficiently large lattices at low temperatures can be simu- 
lated. 



III. DIAGONAL DISORDER 



to a finite DOS at the Fermi level and to a broadening 
of the DOS. For large A the DOS is dominated by the 
disorder spectrum, giving a bandwidth proportional to 
A, and due to normalization a value at the Fermi en- 
ergy proportional to 1/A. This relation is also observed 
for the compressibility at U — at a finite temperature 
(Fig. 2). 

N = 4x4 (3 = 3 U = 4- 




Diagonal disorder describes the idealized situation of a 
random alloy with negligible lattice distortions but vary- 
ing values of the chemical potentials of the constituents. 
While the repulsive interaction U tends to induce singly 
occupied sites, i.e. local magnetic moments, a wide spec- 
trum of random potentials has the opposite effect because 
electrons tend to doubly occupy the lower potentials. In- 
tuitively one would expect that a disorder strength A of 
the order of U may be sufficient to destroy AFLRO . This 
has actually been observed in the limit D — > oo Q where 
the disorder effects are exactly treated by the "coherent 
potential approximation" . |l6| Since the magnetic mo- 
ment formation is a local effect we expect qualitatively 
the same behavior in D = 2. Fig. 1 shows the local 
spin-spin correlation fuction, C(0, 0), on a 4 x 4 lattice 
at U = 4 as a function of A. C(0,0) decreases monoton- 
ically with A and reaches the non-interacting value, 0.5, 
about A = 2U. This local effect is indeed independent 
on dimensionality as it is also observed in D = 1 |i"o| and 
D = oo. [|| Similar behavior is seen in the spin-spin corre- 
lation function at the largest separation, C(2, 2), and the 
AF structure factor S(tt,tt)/N. They are slightly more 
stable than C(0,0) for small disorder and decrease more 
rapidly for A > 4, also reaching their non-interacting 
values of and 0.0508, respectively about A = 21] . 

More interesting is the behavior of the compressibility 
k (Fig. 2). In a Fermi liquid in the limit T — > 0, n is 
equal to the one-particle density of states (DOS) at the 
Fermi energy. For the pure tight binding model without 
interaction and disorder, n diverges logarithmically with 
T at half-filling due to perfect nesting. |17| Turning on 
the disorder removes the van-Hove singularity, leading 



FIG. 1. Local moment C(0, 0) [triangles], longest range 
spin correlation C(2, 2) [squares] and antiferromagnetic struc- 
ture factor S(tv,tt)/N [circles] on 4x4 lattices at U = 4 and 
T = 1/8, as a function of site disorder A. The noninteracting 
values of each quantity have been subtracted out. 



N = 4x4 (3- 




FIG. 2. Compressibility k at U = [triangles] and U — 4 
[squares] on 4x4 lattices and T = 1/8, as a function of site 
disorder A. The very different behavior at A = reflects the 
suppression of the divergent density of states in the noninter- 
acting limit and the opening of a charge gap. 
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On the other hand, when the interaction strength U 
is nonzero, for small values of A the compressibility van- 
ishes due to the charge gap induced by antiferromagnetic 
ordering. AFLRO is strictly present only in the ground 
state, but for a finite lattice the AF correlation length at 
a finite temperature can exceed the lattice size. [fl3|| Even 
with a fully established charge gap, k is always finite 
at finite temperature, however exponentially suppressed 
(Fig. 2). k starts to increase significantly for A > 2 and 
reaches a maximum at about A = 6. Since the k vs. A 
curve for U — 4 approaches the non-interacting curve 
for large A the finite compressibility is quite likely not 
thermally activated but due to the closing of the charge 
gap. Further, k becomes finite for relatively small val- 
ues of A < 4 where the AF correlations are still strong. 
Although the lattice size is much too small for a definite 
conclusion this indicates that the charge gap may close 
within the AF phase. Evidence for a disorder induced AF 
metal at half filling has previously been found in D = oo 
and also in a Hartree-Fock treatment in D = 3. |7]] 
Note, however, that in the presence of disorder a finite 
compressible does not necessarily imply metallicity be- 
cause of localization effects. 

Fig. 3 shows why we cannot study large lattices and 
low temperatures in the case of diagonal disorder. The 
averaged values of the sign of the product of determi- 
nants, ((cr)), vanishes rapidly with spatial lattice size N 
and inverse temperature (3. A value ((a)) smaller than 
about 0.2 precludes reliable simulations due to a vanish- 
ing signal to noise ratio in the data. 

U = 4- A=6 
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FIG. 3. Average sign ((a)) for the site disordered problem 
at 17 = 4 and A = 3 as a function of inverse temperature (5. 
Circles are 4x4 lattices, while squares and triangles are 6x6 
and 8x8 lattices respectively. This sign problem is absent in 
the case of hopping disorder at half-filling. 



IV. OFF-DIAGONAL DISORDER 

In the case of random hopping elements restricted 
to near neighbor sites on a bipartite lattice there is no 
"minus sign problem" at half filling. Therefore we can 
do a much more detailed analysis of the phase diagram. 
We study square lattices with periodic boundary condi- 
tions up to the size N = 100 (N x = 10). For a given 
disorder configuration, 500-700 Monte Carlo sweeps were 
performed for equilibration followed by 1000-1500 mea- 
surement sweeps. Then, all measured quantities were av- 
eraged over 10-20 different disorder configurations. The 
disorder average is the main source of the statistical er- 
rors. 

As in the previous section, we first study the spin-spin 
correlations as a function of disorder strength A for a 
given lattice size. Fig. 4a shows C(l) over a path in real 
space for different values of A on a 6x6 lattice. AF cor- 
relations are present for all values of A. They are only 
slightly reduced for A < 0.8 and much more significantly 
reduced when A ~ 1.6. Larger lattices show a quite sim- 
ilar picture (Fig. 4b). The local moments v/ C(0, 0) are 
apparently stable for off-diagonal disorder, as is shown 
in Fig. 5 where the behavior of C(l) for 1 = (0, 0) and 
(N x /2, N x /2) and the AF structure factor S(tt, n)/N are 
plotted as a function of A. Unlike the case of random 
site energies (Fig. 1), C(0, 0) is almost unchanged by A. 
However, measures of the long range order are strongly 
affected. While the longer range correlations break down 
for strong disorder (A ~ 2) they slightly increase for 
small A. This slight increase between A = 0.0 and 0.2 
might be due to an enhanced averaged Heisenberg ex- 
change coupling 

(T..)- 4{ ® -J + A2 fl2l 
(Jij)- — -Jo + ^j, (12) 

with Jo = 4:t 2 /U for the not disordered case. However 
this is essentially a strong coupling argument (Jo <C i) 
while the present value is only Jo = t. 

In order to answer the question if there is AFLRO 
in the ground state we calculate C(N x /2,N x /2) and 
S(tt, 7r) for different lattice sizes and extrapolate the re- 
sults assuming a finite size scaling according to (111]). 
For the present lattice sizes (up to 10 x 10 sites) the 
magnetic correlations are saturated at a temperature of 
about T = 1/10 where the finite system is essentially in 
its ground state. For disorder strength A < 1.2 both 
C(N x /2, N x /2) and S(tv, 7t) extrapolate to a nonzero or- 
der parameter M at N x — > oo (Fig. 6). The values for 
M are obtained by a least-square fit of the data with 
N x > 6. For N x = 4 there are apparently deviations 
from the scaling (11). Within the statistical error the 
independent extrapolations lead to the same value of M. 
For A = 1.6, however, there is no long range order. 
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FIG. 4. Spin-spin correlations as a function of separation 
on (a) 6x6 and (b) 10x10 lattices at T = 1/10 and U = 4. 
The value shown at (0, 0) is (7(0, 0) - 1/2. 

The values of M as a function of A are shown in Fig. 7 
where the value for A = is taken from the literature 
|p3| . M is apparently stable for small disorder strength 
A < 0.8 and then decreases and eventually vanishes 
about A « 1.4. The slight increase of magnetic corre- 
lation for small disorder at a given system size (Fig. 5) is 
not observed in the values of the order parameter. 



FIG. 5. Local moment C(0, 0) [circles], longest range spin 
correlation C(5, 5) [squares] and antiferromagnetic structure 
factor S(7r,7r)/N [triangles] on 10x10 lattices at U — 4 and 
T — 1/10, as a function of site disorder A. The noninteracting 
values of each quantity have been subtracted out. 

The physical reason for the destruction of AFLRO is 
not obvious for this type of disorder. There is no mag- 
netic frustration and also no destruction of local mo- 
ments. One possible approach to understanding the tran- 
sition is a weak-coupling analysis. Within the RPA the 
magnetic susceptibility of the interacting system xrpa 
is expressed in terms of the noninteracting value \o : 



Xj?.pa(tt,7I") 



XoO,7r) 



1 - Uxo(n, 7r) 



(13) 



The perfect nesting relation of the pure tight binding 
model without interactions and disorder: e(q+ (ir, 7r))) = 
— e(q) at half filling leads to a divergent susceptibility 



Xo(7r,7r) oc (InT) 7 . 



(14) 



For the square lattice at half filling we have 7 = 2 due to 
the logarithmic van-Hove singularity in the DOS at the 
Fermi energy. In general, for bipartite lattices without 
this DOS singularity like the simple cubic lattice the ex- 
ponent becomes 7 = 1. Since xo diverges at T = the 
RPA predicts that the system will order at any finite U. 

One effect of the disorder is removing the perfect nest- 
ing relation of the pure tight binding model, reducing 
Xo and leading to the destruction of AFLRO within this 
RPA analysis. However, off-diagonal disorder does not 
remove the bipartite structure of the lattice, so we ex- 
pect Eq. ( jl4| ) still to hold with 7 = 1 although q is no 
longer a good quantum number. Since there exists no 
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FIG. 6. Finite size scaling analysis of the AF structure fac- 
tor and long range spin correlations. These quantities, appro- 
priately scaled, are plotted as a function of the inverse linear 
lattice size 1/N X . A non-zero extrapolation to 1/N X = indi- 
cates AFLRO. The extrapolated values should be identical. 
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FIG. 7. The order parameter M inferred from Fig. 6 as 
a function of A. For A < A c » 1.4 there is AFLRO. For 
A > A c « 1.4 there is a transition into a disordered phase. 

analytic solution for (x(q)) for the disordered lattice even 
without interaction, we calculate it numerically for finite 
lattices. xq( 7T i 7T ) is indeed well described by a logarith- 
mic fit (Fig. 8). In fact, we further observe that in the 
disorder-averaged RPA susceptibility XRPA (""j 7r) f° r fi- 
nite U one can to a very good approximation replace 
x(7r, 7r) for a given disorder configuration by its average 
value (x(tt, tt)): 



Xapprox. (jt) ft) 



X(7T,7r) 



1 - ^X(7r,7r) 

(X(7T,7T)) 

l-C/{ X (7r,7r)) 



(15) 



In particular the estimated Necl temperatures for given 
U obtained from the divergence of xrpa(tt, tt) and 
,(7r,7r) are identical within the errors due to dis- 
order averaging (Fig. 8). We observe that xrpa(<i) for 
different momenta q ^ (tt, tt) is always more strongly 
suppressed by disorder than xrpa{tt, tt). We conclude 
that off-diagonal disorder does not remove the low tem- 
perature divergence of the AF susceptibility. Thus an 
explanation of the transition based, for example, on the 
inducement of a finite U cr it > 4 as the disorder is in- 
creased is not viable. 

While the above considerations might be suited in the 
case of weak interaction, we will now consider the oppo- 
site case of U ^S> t. Here, the half filled Hubbard model 
with random £y at half filling becomes equivalent to the 
disordered, but unjrustrated, AF spin-1/2 Heisenberg 
model with random couplings Jy, see (|l^). This model 
has been studied in D = 2 by a Quantum Monte Carlo 
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FIG. 8. The behavior of the noninteracting susceptibility 
Xo (circles) with temperature T in the presence of disorder. 
Xo still diverges (logarithmically) as T — > 0. The full line is a 
fit with xo = 0.054 — 0.245 In T. Also shown are the inverse 
RPA susceptibility XfiPA( 7r i 7r ) (triangles) and the approxi- 
mation Xapprox. i 77 j 71 ) (squares). They are indistinguishable 
within the statistical errors and vanish at the same critical 
temperature. The dotted line is a guide to the eye only. 

technique in the case of a binary distribution of Jy G 
{J\,J2~\ II- This numerical work incorporated disorder 
in two ways- both through the difference of the exchange 
energies A = | Jy — J%\ and also through different con- 
centrations of strong and weak bonds. Depending on the 
concentration p of, say the J\ bonds, and the difference 
A, an AF ordered and a disordered phase were identi- 
fied. The basic physics is that spin singlets form on the 
strong bonds, driving the formation of a paramagnetic 



phase, as has been discussed by Bhatt and Lee. p9 
We have made a quantitative comparison of our data for 
the Hubbard model at U — 4 with this strong coupling 
theory which incorporates spin degrees of freedom only. 
We first note that the result of ||], for the phase bound- 
ary in the p — A plane is reasonably well described by 
setting the variance of J 1} , 



P(1-P)A 2 



(16) 



to a critical value v — v c . |2lJ That is, the fact that there 
are two distinct types of randomness appears irrelevant- 
they can be modeled together in a simple way by their ef- 
fect on the variance of Jy . The fit to the calculated phase 
boundary is best for v c w 0.40 to 0.42. The variance of 
J in the Hubbard model with disorder distribution (JsJ) is 
independent of U and reads: 



7 



2 \2 



60i 2 A 2 



A 4 + 40i 2 A 2 + 80 



(17) 



u has its maximum of 5/9 at A = 4yfEt. If we now apply 
the same criterion, v = v c , as in the Heisenberg model 
we obtain a critical disorder strength of A c = 1.7t to 
1.8t. This estimation is surprisingly close to the critical 
A obtained by finite size scaling at U = 4 (Fig. 7) in spite 
of the fact that U = 4, corresponding to Jo — 1 is n °t 
in the strong coupling limit of the Hubbard model. fl§| l 
The agreement might be due to the fact that already for 
U = 4 the double occupancies are strongly suppressed 
and the density of local moments has reached about 75% 
of the maximal value 1.0 independent of disorder strength 
(see Fig. 5). Further, A c is small compared to 3J7 so that 
the averaged (Jy) is onl y s lightly enhanced (by a factor 
of ~ 1.25) according to ( fl2| ) which is a consistency check 
of the arguments above. 

To characterize the properties of the disordered phase 
(A > 1.4) we study the temperature dependence of the 
uniform spin susceptibility %(0, 0) = f3S(0, 0) (||) and the 
charge compressibility k ( |l0| ) . In the non-interacting case 
both quantities are identical, and diverge logarithmically 
in D = 2 for A = 0, as mentioned before. flTj This sin- 
gularity is removed by interactions, and x(0j 0) reaches 
a maximum at a finite temperature, approaching a fi- 
nite value at T = 0. [|l3| Disorder suppresses the non- 
interacting susceptibility too, but for a different reason 
because it removes the van-Hove singularity in the den- 
sity of states (Fig. 9a). In the interacting case however, 
disorder has the opposite effect of enhancing x(0, 0) if the 
disorder is strong enough (A = 1.6). For A < 0.8 there is 
no significant effect on x(0, 0) due to disorder at the low- 
est temperature under consideration (T = t/10). In the 
present temperature regime at A = 1.6, x(0, 0) decreases 
monotonically with T, a behavior also found in the dis- 
ordered phase of the random bond Heisenberg model || . 
Stronger disorder (A = 2.4, see Fig. 9b), where some of 
the hopping elements become negative, leads to further 
enhancement of x(0, 0). In systems with longer range 
interactions a low temperature divergence of x(0, 0) has 
been predicted. |19| This divergence is explained by the 
presence of localized moments, i.e. spins on lattice sites 
which arc effectively decoupled from the rest of the sys- 
tem. Such isolated sites are particulary likely in the case 
of a special type of correlated hopping randomness with 
Uj = %iX$ where the site variables {xi} are randomly 
distributed. j2^] In contrast to the uncorrelated distribu- 
tion (^3j), here a site i with a small variable Xi is weakly 
connected to all its neighbors. We consider a bimodal, 
symmetric distribution of site variables: 



P( Xi ) 



A' 

T 



6(x\ 



(18) 



with A' < It. This distribution corresponds to a binary 
alloy where the hopping amplitudes t» depend only on 



the (three possible) combinations of atoms on neighbor- 
ing sites, ijj are always positive and can take the values 
t - A'/2, y/t 2 - (A'/2) 2 , and t + A'/2. The average 
Heisenberg coupling remains (Jy) = J = At 2 /U inde- 
pendent of A'. We find that for A' = 1.6, x(0,0) 



LO 

CD 




0.6 



T 




0.6 



T 



FIG. 9. Uniform susceptibility x(0, 0) vs. T on a 8 x 8 lat- 
tice, (a) x(0,0) for U = 4 and A = 1.6 (triangles); x(0,0) 
for U = 4, T = t/10, and A = 0.0, 0.4, and 0.8 (circle) 
are indistinguishable within the error-bars. Also shown is the 
non-interacting n = x(0, 0) for A = 0.0 (dotted line) and 
A = 1.6 (full line). The compressibility k for U — 4 and 
A = 1.6 (squares) vanishes exponentially at low temeratures. 
The key point of this data is that while randomness suppresses 
the susceptibility at U — 0, it enhances it for U nonzero, (b) 
k = x(0,0) for U = A = (dotted line) and x(0,0) for 
U — 4, A = 1.6 (triangles) are again shown as in (a). How- 
ever, here they are now compared with x(0, 0) for A = 2.4 
(squares) and for the case of correlated hopping randomness 
with A' = 1.6 (circles) (see text). 
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increases for lower temperatures and is enhanced by al- 
most one order of magnitude at T = t/10 compared to 
the case without disorder. Although the statistical er- 
rors and limited temperatures do not allow to extract 
the functional dependence on T, this dramatic increase 
supports the existence of a finite density of localized mo- 
ments. 

As in the interacting case without disorder, the com- 
pressibility k, shows for A = 1.6 activated behavior with 
temperature (Fig. 9a) with an estimated charge gap of 0.5 
p3| . A quite similar behavior is found even for stronger 
disorder A = 2.4 and also for the correlated hopping 
randomness with A' = 1.6 (not shown): n is very small 
(~ O(10~ 3 )) and decreases strongly for lower T consis- 
tent with an activated behavior. Hence the systems re- 
mains uncompressible even in the absence of long range 
magnetic order. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have presented Quantum Monte Carlo 
calculations for the magnetic phase diagram of the dis- 
ordered repulsive Hubbard model at half-filling. When 
the randomness is in the site energy, the sign problem 
precludes the study of large systems at low temperatures 
necessary for a proper finite size scaling analysis. How- 
ever, we are able to see on 4x4 lattices the destruction of 
local moments and longer range spin-spin correlations. 
We also studied the compressibility function of 

disorder strength. The incompressible state with n = 
at nonzero U in the absence of disorder is replaced by a 
finite k at A « 2, indicating a disorder induced closing 
of the charge gap. 

When the randomness is in the hoppings iy a particle- 
hole transformation proves that the product of determi- 
nants which gives the Boltzmann weight is positive at 
half- filling. That is, there is no sign problem at arbi- 
trarily low temperatures and large lattices. In this case 
we were able to do a finite size scaling analysis of our 
results and determined A c » 1.4 for the destruction of 
antiferromagnetism. This value was in reasonable agree- 
ment with that obtained for the disordered Heisenberg 
model ||. The paramagnetic phase is driven by singlet 
formation. Unlike the random site energy case, the lo- 
cal moments remain well-formed through the transition. 
The strong enhancement of the uniform susceptibility at 
low temperatures, in particular for the correlated hop- 
ping randomness, indicates the formation of effectively 
localized magnetic moments. The compressibilty is ap- 
parently activated even for strong disorder where there 
are no long range magnetic correlations. 

We are currently studying the dependence of the den- 
sity of states in this model as it has recently been sug- 
gested that this may provide an alternate way to describe 
transitions in the Anderson-Hubbard model and its ex- 
perimental realizations. p3| 
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